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Abstract 

We consider the geometric entropy of free nonrelativistic fermions in two di- 
mensions and show that it is ultraviolet finite for finite fermi energies, but 
divergent in the infrared. In terms of the corresponding collective field the- 
ory this is a nonperturbative effect and is related to the soft behaviour of 
the usual thermodynamic entropy at high temperatures. We then show that 
thermodynamic entropy of the singlet sector of the one dimensional matrix 
model at high temperatures is governed by nonperturbative effects of the 
underlying string theory. In the high temperature limit the "exact" expres- 
sion for the entropy is regular but leads to a negative specific heat, thus 
implying an instability. We speculate that in a properly defined two dimen- 
sional string theory, the thermodynamic entropy could approach a constant 
at high temperatures and lead to a geometric entropy which is finite in the 
ultraviolet. 



Recently the entropy of entanglement between different regions of space 
in quantum field theories have been intensively studied The moti- 

vation for this is its direct connection to the question of information loss 
due to black holes and black hole entropy f7fl-112||. A significant feature of 
this entanglement entropy, or "geometric entropy" is that it is ultraviolet 
divergent in typical field theories. This has been interpreted to imply that 
at least at the semiclassical level information loss due to the formation of 
a horizon is inevitable in quantum field theories. The divergence of the en- 
tropy is a reflection of short distance singularities in quantum field theories. 
Alternatively |I3]|-||16|| the divergence is related to the behaviour of the usual 
thermodynamic entropy at high temperatures since, as we shall see below, 
the geometric entropy effectively involves an integral of the thermodynamic 
entropy density over all temperatures. 

One may hope that in string theories this divergence disppears because 
of a soft ultraviolet behaviour JXX[] . However, to leading order in the string 
perturbation expansion, the thermodynamic free energy of a string is equal 
to the sum of the free energies of the physical modes of the string and one 
would obtain the same divergence in each term of the sum. Furthermore, 
unlike in a field theory of a finite number of fields, the thermodynamic free 
energy of free strings is itself divergent at the Hagedorn temperature. This 
is an infrared diveregence and might signal an instability of the theory and 
one might obtain a finite answer once one takes into account interactions and 
shift to a stable vacuum [lTj. As argued in [13, 14| the geometric entropy in 
string theory is afflicted by this Hagedorn transition and what appears as an 
ultraviolet divergence in each term of the sum over all string modes may be 
interpreted as an infrared problem in the full answer. 

The need to include string interactions calls for a formulation of the prob- 
lem in some well defined and tractable string field theory. While this appears 
to be an almost impossible task at present, there is one string theory where 
a tractable nonperturbative formulation exists, at least for some bulk quan- 
tities. This is the two dimensional string defined via the one dimensional 
matrix model [IS]. The singlet sector of the matrix model [] may be written 



as a two dimensional collective field theory of the density variable. The fluc- 
tuations of the collective field represent a massless particle which is the only 



1 As we will see soon the nonsinglet contributions are irrelevant for a calculation of the 
geometric entropy in the ground state 
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propagating degree of freedom of the two dimensional string Q and the cou- 
pling is proportional to the inverse of the fermi energy. For nonperturbative 
considerations it is better to write the model as a field theory of nonrelativis- 
tic fermions in the presence of an inverted harmonic oscillator potential and 
no self-interactions. The idea, then, is to consider a geometric entropy in 
this model and use exact nonperturbative answers to understand stringy ef- 
fects. Hopefully this will teach us something about higher dimensional string 
theories as well. 

In this paper we take the first step in this program. We first consider the 
problem of free nonrelativistic fermions in two dimensions. By constructing 
the explicit expression for the ground state wave functional and the corre- 
sponding geometric density matrix we argue that the geometric entropy has 
no ultraviolet divergence for finite fermi energies. In terms of the collec- 
tive field theory this is a nonperturbative phenomenon. In fact the result 
follows from the softer behaviour of the ordinary thermodynamic entropy 
at high temperatures (compared to relativistic fermions). We then consider 
the the high temperature limit of the "exact" expression for the thermody- 
namic partition function of the singlet sector of the one dimensional matrix 



model [21|. We will show that the genus expansion breaks down at high 
enough temperatures and hence the geometric entropy of this string theory 
is essentially nonperturbative. A naive high temperature expansion leads 
to a regular behaviour of the entropy, but the specific heat turns out to be 
negative, signifying a nonperturbative instability. 

Consider a theory of fields 4>{x,t) in 1 + 1 dimensions. The ground state 
wave functional may be written as ^/ [4>l, <Pr], where (pL (0#) denotes the 
field (f) for x < (x > 0). The density matrix which gives expectation values 
of operators localized in the x > region is 

P (<f) R , <f>' R ) = J V(P L ^o[(pL, M ^ [4>L, <f>' R ] (1) 

As shown in |Q, Q the quantity Trp n may be represented as an euclidean path 
integral over a cone with a deficit angle 2n(l — n). The geometric entropy 

2 The relationship between the collective field and the massless scalar of the effective 
field theory is rather subtle and not completely clear at this moment. See |23|, ||| for a 
recent discussion. 
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may be then written using a "replica trick" as 

S g = -plog p = [(1 - ™^)logTr p"] n=1 (2) 

For relativistic systems this establishes the equality of the geometric entropy 
with the usual thermodynamic entropy of the field in Rindler space at the 
Rindler temperature, and hence the quantum correction to the entropy of a 
large mass black hole. For nonrelativistic systems the Rindler hamiltonian 
would depend on the Rindler time. However the geometric entropy defined 
above makes sense and Trp ra is still a path integral on a cone. 

Consider a system of N F free nonrelativistic fermions contained in a box 
: — L < x < L. The dispersion relation for the single particle states is given 
by e(k) = \k 2 where k is the (spatial) momentum which is quantized as 
k — with integer or zero n. However, below we will often use continuum 
notation. The second quantized fermion field operator can be expanded in 
terms of quasiparticle operators b(k) and b + {k) 

POO Ah 

#M) = / ^e- ikx [b(k)9(\k\-k F )+b + (-k)9(k F -\k\)} (3) 

and similarly for The operators b, b + satisfy the standard anticommu- 
tation relations {b(k, t), b + (k', t)} = S(k — k'). The ground state is then 
described by the filled Fermi sea b(k)\0 >= 

In terms of the coherent states of the 6-oscillators 

b(k)\b(k) >= b(k)\b(k) > <b{k)\b + {k) =b(k) <b(k)\ (4) 

where b(k),b(k) are grassman numbers, the ground state wave functional is 
given by 

*„[&(*), Kk)\ = exp[-- j ^ -b{k)b{k)\ (5) 

The grassmann fields appearing in the path integral are, however, not b(k) 
and b(k). Rather, they are the combinations 

ijj(k) = b(k)e(\k\ - k F ) + b(-k)6{k F - \k\) (6) 

and similarly for ip(k). The fourier transforms of these fields are the original 
fields ip(x). For reasons which will be clear in a moment we will use a new 
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field x(q) — + q). Then the wavefunctional (]5|) becomes 

*o = exp r ^x(ff)x(ff) - 2 /° P-X(q)x(q))\ (7) 

The inverse fourier transform of x(q)i denoted by x{ x ) is related to the 
original field ip(x) by x( x ) = e lkFX ip(x). Since this relationship is local we 
can compute the geometric entropy using these fields. 

In (0) the first term may be written as a integral over position space of a 
local quantity x( x )x( x ) an d do not contribute to the geometric entropy The 
only momentum modes which contribute to the geometric entropy are those 
in the filled Fermi sea as in the second term in (|^). The fermi momentum 
thus roughly acts as an ultraviolet cutoff and this is the essential reason why 
the geometric entropy turns out to be finite. 

We can proceed similarly for a relativistic Weyl fermion. For example 
for right moving fermions the fermi momentum is zero and the Dirac sea 
consists of all negative momenta. It is easy to see that the wavefunctional is 
the expression (|7]) where the limit of integration —2kp in the second term is 
replaced by — oo. Clearly, in the limit of large kp the nonrelativistic expres- 
sion reduces to the relativistic expression. This simply reflects the fact that 
the excitations very close to the fermi level behave as massless relativistic 
particles with the velocity of light replaced by the fermi velocity. 

To compute the geometric entropy we now introduce the following tech- 
nique which may be trivially generalized to other situations (e.g. relativistic 
bosons, fermions etc.). We first expand the field eigenvalues in terms of modes 
which are localized to the left region (x < 0), /z,(w) and those localized to 
the right region (x > 0), /r(w) as follows 

j fOOnf /"OO rp 

X (x) = -f=[0(x) / duo{-r^f R {uo) + 6{-x) / du(—)-*"f L (u,)] (8) 

\ \x\ 00 ^ 00 ^ 

where a denotes the lattice spacing (or some other ultraviolet cutoff). The 
factor of — 7= arises from the dimension of the field y(x) under rescalings of 

x. To obtain the expression for x + ( x ) replace fi, /r by fi, and u by — uj 
in the integrand of (Q). 

The fields x(q) an d X + (<?) may be now expressed in terms of fp and Jr. 
We get the following expressions, all valid for q > 0. 

/OO rjijj 
— [ie™f L (u) + f R (u)]G(q,u) 
-oo Z7T 
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/°° rlii) - - 
— [f L (oj)+ie-™f R (u)]G(q,-u) 
-oo Z7T 

Xi-q) = [°° ^[f L (u;) + ie™f R (u)}G(q,u) 

Xi-q) = p-[ie^h(u>) + f R (uj)]G{q,-u>) (9) 



/•OO (inn q. 1 ?, 1 . 

G(q,u)= ^e iqx (-y w = a^{—Y^T{--iu) g>0 (10) 
Jo \/x a qa 2 



where we have defined 

dx ..„., ./•.,_, . i i. 1 

Note that G(— g, w) = i e nul G(q,u), which has been used in writing 

The ground state wave functional may be now rewritten explicitly in 
terms of f L ,f R : 

*o[/l, fn] = exp [--(/(oo) + /(oo) - 21(2**))] (H) 



J (7) = / 7rx(q)x(q); = I 7r-x(q)x(q) = I 7rx(-q)x(-q) 



where we have defined 

[^mx{^ j Jo 2n 

(12) 

In evaluating 1(7) we need to perform the integral 

F 7 {u>, co') = £ ^G(q, u) G(q, -to') (13) 

In the presence of a finite box size 2L the lower limit of the q integral is really 
Using (Q) one may write the wave functional explicitly in terms of fi and 
f R and obtain the density matrix for fields localized in the x > region by 
simply functionally integrating over f L , f L 

p(f R , f R ) = J vj L vh MfR, h] Mfk, h\ (14) 

It is easy to evaluate the integral F 7 (u;,a/) for 7 = 00. The result is 

Foo(uj,ui') = 2ni sech(7ru;)5(u; — to') (15) 
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From this we can make the first consistency check on our technique. Using 
( |15D one gets the wave functional for relativistic fermions 

% el = exp [-2-Kfdu [tanh Truj(f R (u)f R (uj) - f L (w)f L (w)) 

+i sech ™(f L {u)f R (u) - f R {u)f L (u))] (16) 

After a trivial rescaling of the fields by a constant | this is exactly the answer 
obtained by euclidean path integral methods in |5[]. As discussed above, (|IT)|) 
is also the answer for nonrelativistic fermions exactly at hp = oo. 

The density matrix and its various powers in this relativistic limit may 
be obtained as in || and the geometric entropy may be calculated using the 
replica trick to obtain the result 

S g ~ log(L/a) (17) 

This expression has an ultraviolet as well as an infrared divergence. This 
result may be understood in a very simple way. The geometric entropy is a 
measure of the entanglement of the modes fp and f R . It is clear from the 
wave functional that the modes /l(^) and f R (u) mix for all values of uj. 
This means that all the cu-modes contribute to the geometric entropy, which 
should be then proportional to the total number of uj's. From (|8]) it follows 
that for large L the number of allowed cj's is proportional to log(L/a) - hence 
the answer (|P7D. 

For small values of kp the integral l{2kp) may be approximated as 

I(2k F ) ~ (2k F - 1) x(0) X (0) (18) 

It may be easily seen that for large N = — the function G(0, ui) is peaked 
at oj = 0, falls to zero at uj ~ and then oscillates rapidly around zero. 
This means that the mode x(0) expressed as an integral over uj essentially 
receives contribution from a few values of to around uj = 0. To the lowest 
order we can then approximate the wave functional by 

/oo _ _ 

duj{f L {uj)f L {uj) + f R {u)f R (u)) 
-oo 

+21og (^) [-f R (0)f R (0) - /i(0)/i(0) 
+if R (0)f L (0) - i7l(0)/x(0)] (19) 

3 It may be easily verified from the definitions that a rescaling of fields by a constant 
does not change the geometric entropy, as one would expect. 
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In ([19]) the original integrals over u have been replaced by sums and the 
resulting factors of logiV have been suitably absorbed by a redefinition of the 
fields. 

In ( |19"D the left and right modes mix only for u = 0. For small but finite 
kp only a few modes around uo = mix, and these alone contribute to the 
geometric entropy. Let us evaluate the entropy due to the to — mode alone. 
The final result will be this contribution multiplied by a factor of order unity. 
The unnormalized density matrix for this mode is easily seen to be 

pUrJ'r) =exp[- (l-~sech 2 V )(f R (0)f R (0) + f' R (0)f' R (0)) 

+ isech 2 r / (/ /? (0)^(0) + &(0)/«(0))] (20) 



where we have defined 



log^ 

sech 7] = - — (21) 

' log^ + 2 v ' 



By a redefinition of the fields f R — > (2tanh 77)2 f R one may rewrite the 
density matrix in a form which facilitates the calculation of Trp n for any n. 
The final result for the geometric entropy is 

S g = 21og(2cosh 7]) - 2r] (tanh rj) (22) 

For kpL ~ 7r one has i] — > ±oo and the leading order expression for the 
geometric entropy follows from (p2|) 



^ 0) = (log^) 2 [l + log 2 - 2 log (log^)] (23) 

7T 7T 

The total geometric entropy for small kp is multiplied by a factor of 
order unity. As advertized above the result does not involve the ultraviolet 
cutoff essentially because only a few modes contribute to the entropy, rather 
than all of the log iV modes. In a similar fashion we expect that for finite kp 
the role of the ultraviolet cutoff is replaced by kp. 

There is an alternative way to view the quantity Trp n . Consider dividing 
up the region x > into small cells. In the thermodynamic limit the path 
integral which represents Trp n would be a product of path integrals for each 
individual cell. However the path integral for the cell centered at the point x 
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is the standard thermodynamic partition function at a temperature T{x) = 
Thus the geometric entropy may be obtained by simply calculating the 
thermodynamic entropy density at a temperature T(x) and integrating the 
result from x = to x = oo. This is the procedure used to compute the 
genus one contribution to the entropy of strings in |L3|, |i"4| . 

The ultraviolet divergence of the geometric entropy in relativistic field 
theories may be now understood in terms of the behaviour of the standard 
thermodynamic entropy at high and low temperatures. For example, a free 
massless boson in d space dimensions at temperature T = \ has an entropy 
density s = -A. This diverges for low /3, or high temperature for all d. In 
our problem we have to put (3 = 2ttx and integrate this entropy density over 
x, so that there is a diveregence of the geometric entropy from the lower 
limit of integration. This is the ultraviolet divergence in the entropy. In 
the corresponding Rindler problem, the divergence arises because the local 
temperature is very high near the horizon and the corresponding contribution 
to the entropy density is large. For d — 1 there is an additional infrared 
divergence coming from large x. Thus using an ultraviolet cutoff e and an 
infrared cutoff L in the integral over x one has S ~ -j=t for d > 1 and 
S ~ log© for d=l. 

Consider now our system of Np free nonrelativistic fermions in one spatial 
dimension in the grand canonical ensemble with an inverse temperature (3 = 
^ and fugacity z. The fugacity is determined in terms of Np by Np = 
zd z log Z (where Z is the partition function) which leads to 

N F = Uj^-h (z) h (z) = zdjs (z) (24) 
where the function f 3 is defined as 

J 2 

1 f°° 9 

fs(z) = — dx\og(l + ze- x ) (25) 

2 J-oo 



The expression for the entropy follows from usual thermodynamics and is 
given by 

3f 3 (z) 

The function fi (z) is a monotonically increasing function of z. Thus for 
large values of (i.e. low temperatures or high densities) one may use 
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the standard expansions of the functions f r {z) for large z [HJ to obtain the 
leading order term in the entropy density 

Np 7T 2 13tt 4 tt 8tt 3 

S ~ lp61og z + 360(log zf + " ' J ~ 6/?v^ + 45 /3 3 fc| + " ' ^ ^ 

The leading term is the same expression for a free relativitsic boson or a 
free relativistic fermion in two dimensions, apart from a factor of y/2ep. The 
reason is simple. For a given temperature and size of the system, large 

values of mean large fermi momenta. In this case the excitations are 

restricted to particle-hole excitations near the fermi level. The dispersion for 
these excitations are given precisely by E(k) = y/2ep\k\ which is the same as 
that of a massless relativistic boson in two dimensions. The factor of y/2ep 
in the dispersion relation explains the same factor in the entropy density. 
What is more significant for our considerations is the fact that the same 
result applies for a given \ip and size, but low temperatures. 

On the other hand, for small values of - one has to use the power 



series expansions of f r (z) for small z [15) and one easily gets for the entropy 
density 

s = ~f ~ M— ^) + 0{F) (28) 

To obtain the geometric entropy one has to integrate s(2irx) over all positive 
x. If we use the low temperature expansion we would get an ultraviolet 
divergent answer from the behaviour of the integrand near x — 0. In fact it 
may be easily checked that the lowest order answer agrees completely with the 
direct calculation of the geometric entropy of relativistic fermions discussed 
above. However to treat the region x = we have to use the high temperature 
behaviour in (p5[) . This shows that the integrand has an integrable singularity 
at x = whereas there is a logarithmic divergence coming from large values 
of x. Thus the geometric entropy is finite in the ultraviolet. This is in marked 
contrast with relativistic bosons or fermions. Essentially as one approaches 
the point x = the temperature becomes large compared to the scale set by 
the fermi energy. Thus the fermi energy provides a cutoff to the relativistic 
behaviour for large values of x. 

The true significance of this result can be appreciated if one rewrites the 
model in terms of the collective field theory of the density of fermions p(x, t). 
Consider a more general system of iVV fermions in an external static potential 
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V(x). The density p(x,t) may be expanded around the "classical" average 
value 

1 



p(x,t) = po(x) + d x £(x,t) p = -J2(p F -V(x)) (29) 



7T 



where p F is the fermi level. Introduce the "time of flight" variable t(x) = 
\ p^)' Then the hamiltonian which governs the fluctuations £(r, t) is 
given by 

H coll =:\j rfr[n|+(9 T o 2 -^(n,(9 T on,+^(9 T o 3 )-4(||-^)] : 

l J ix 2 pg 6 7r 2 op zp 

(30) 

where 11^ denotes the field momentum. Note that the hamiltonian already 



comes in the normal ordered form pOl and leads to finite answers. The 



density thus behaves as a scalar field with in general position dependent 
coupling given by For free fermions V(x) = 0. Then p = y/2p F and 

t = \ The coupling is a constant and equal to j^. F° r Pf large 
compared to typical energies, one has a free relativistic scalar field, which 
is why in this limit we obtained the relativistic answer for the entropy. For 
small hf the theory is strongly coupled and the perturbation expansion does 
not make sense. Thus the finiteness of the entropy demonstrated above is a 
non-perturbative effect in this collective field theory. 

Let us now consider the one dimensional matrix model which provides a 
nonperturbative definition for the two dimensional noncritical string. This 
is defined by the action 

A — X J dm (^[d t M(t)) 2 + V(M(t))] (31) 

where M(t) is a N x N matrix and the potenial function V(M) has to 
be chosen such that it has a quadratic maximum. The detailed form of the 
potential is not important in the double scaling limit. As is well known the 
singlet sector of this model may be written exactly in terms of nonrelativistic 
fermions ip(x,t) where x denotes the space of eignevalues of the matrix M. 
These fermions have no self interactions, but move in an external potential 
V(x). Let us denote the fermi energy by pp. Then as the coupling g = y 
of the theory approaches a critical value g c the fermi level p F appraches p c 
which is the energy at the top of the potential hump. The double scaling 
continuum limit of this problem is given by p = p c — pp ^ and A — > oo 
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with k = A/i = fixed. This is known to describe the two dimensional string. 
In this limit only the quadratic hump of the potential is relevant and one 
has a problem of nonrelativistic fermions in two dimensions in an external 
potential V(x) ~ — x 2 - the inverted harmonic oscillator. The collective field 
theory hamiltonian is given by (j30|) with this potential. This is now a version 
of string field theory and the string coupling is g st ~ -. 

The geometric entropy of this model in the ground state may be calculated 
in principle by the technique discussed above. The expansion for the fermion 
fields are no longer in terms of plane waves, but in terms of the eigenf unctions 
of the Schrodinger operator in the inverted harmonic oscillator potential, 
which are parabolic cylinder functions. The main modification would be to 
replace the plane waves e iqx by parabolic cylinder functions in the definition 
of G(q,uj) in (|K]). The resulting expressions are rather difficult to analyze. 

We have seen, however, that the divergence of the geometric entropy is 
related to the high and low temperature behaviors of the ordinary thermo- 
dynamic entropy. The complete thermodynamic partition function of the 
matrix model is a formidable task since this includes contributions from the 
non-singlet states. However if we are interested in computing the geometric 
entropy for the ground state, which is a singlet, the trace involved in Trp n is 
a trace over singlet states alone. The thermodynamics in the singlet sector 



has been completely solved and we can use the known results of [^1[ 

Let us define A = g c — g and let F denote the free energy. Then the 
chemical potential /z is determined by the equation 

dA I dt _ it (t/n) (irt/nf3) 1 A 

^7 ~ 27 7o T 6 sinh(t/ K ) sinh(7rt/«/3) J + 2^ [ g 2 ~ 7j ( } 



The infinite constant 7 is necessary for (p2|) to reproduce the correct WKB 
expansion in jk but would be unimportant in what follows. The free energy 
F is then obtained by integrating the equation 

r)F 

^ = A 2 (/i-^ c ) (33) 

Recalling that the number of fermions is iV = Ag one may write down the 
expression for the entropy using standard relations in the grand canonical 
ensemble 
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The second term in (K3J) arises because we are considering derivatives with 
fixed \x rather than fixed N. 

The above expressions have an important duality symmetry. From (|3"^) 
it follows that 

dA, n x . dA,n 2 A/3, 

^^^•v 1 (35) 



It then follows from fl3"3] ) that the canonical partition function f3F is invariant 
under this duality transformation. 

The standard genus expansion is obtained by considering k to be large 
and expanding the hyperbolic sine functions in a power series expansion. The 
result for this asymptotic expansion is 

^"^[-logA.+ E^p] (36) 
where the functions f m (/3) are symmetric under the duality transformation 



i 



j| and has the form 



fc=0 

C(m, k) are numbers related to Bernoulli coefficients. The m = 1 in the sum 
in the above expression is the one loop contribution in string theory. The 
corresponding one loop free energy is 

\0g ti p 7T 

F ~mT^ + W (38) 

This result is identical to the free energy of an ideal gas of massless bosons 
in one spatial dimension of size log \x apart from a constant which is in fact 
the jinite one loop correction to the ground state energy of the system. This 
is expected since the only propagating mode of the two dimensional string is 
a massless scalar. The same result is obtained by performing the Polyakov 
path integral in the continuum d = 2 string theory and is in fact the 
result of the following modular invariant integral 



d 2 r ^ P 2 \n — mr 

-^r > exp[ 

t t£ m „ 4ttt 2 



expi- ' jzz ' ] (39) 
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where r is the complex modular parameter on the torus and the integration 
is over the fundamental domain. The term in the sum with m, n = is 
the zero temperature free energy and is separately modular invariant. The 
temperature dependent term evaluates to the second term in ( |38|) . As noted 
this is the contribution from a single massless scalar. This is nevertheless 
modular invariant as it should be since it follows from a string theory Q . 
Note further that unlike in higher dimensions the one loop answer does not 
have a Hagedorn behavior at any finite temperature simply because the 2d 
string has only one propagating degree of freedom. 

The contribution to the geometric entropy from the one loop term alone 
gives the standard logarithmically divergent answer for a free relativistic 
boson in two dimensions. One may regard the entire divergence to be infrared 
since the answer is modular invariant, as argued in [14]]. However the answer 
is really the same as a single boson. The string interpretation then gives a 
relation between the infrared and ultraviolet cutoff in terms of the cutoff in 
the modular parameter Imr. 

Returning to our description of the usual thermodynamics of the ma- 
trix model let us record below the genus expansion for the thermodynamic 
entropy. 

S = -^[lo g/1+ l] + k^ + 9^] + ... < 4 °» 

The asymptotic expansion (^) makes sense for small temperatures, i.e. 
large (3. Indeed, for large (3, |p approaches a constant, and the leading 
correction from all the terms in the sum in fl36|) is of order Jj, as is clear 
from (|37|) . The leading term in the free energy is then seen to be of order jp 
as well and has contributions from all genus. 

We have seen, however, that the ultraviolet divergence of the geometric 
entropy is related to the high temperature behaviour of the free energy. The 
asymptotic expansion of ( |32"D leading to the genus expansion (^) clearly 
breaks down for high temperatures. More specifically for (3 < - one should 
not expand the hyperbolic functions in (|32]) in power series. This means 
that for high temperatures strong coupling effects (in the sense of string 
theory) become important. Indeed the genus expansion in (^) contain ever 
increasing powers of ^ with higher and higher genus giving rise to higher and 

4 See [[161 for related remarks 
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higher powers. It is clearly necessary to look at the high temperature limit 
for arbitrary values of the string coupling g st . 

Before considering the behaviour at high temperatures let us look at 
the behavior at low temperatures,/? >> ~, but independent of the genus 
expansion. This means we expand the factor (nt/ k0) /smh.(nt/ nf3) in fl32| ) in 
a power series, but not the first factor. This gives a power series in for all 
values of k. A straightforward evaluation of the integrals yield the first few 
terms 

f - ~Re*[i(l + ^)] + ^Ba^[i(l + «)] 

'"' +*«)] + ••• (41) 
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where ty(x) denotes the digamma function. It may be checked that if one 
now considers the large k asymptotic expansion of the digamma functions in 
(flip , the result agrees with the large- 6 limit of the genus expansion in fl36l) . 



Since we have duality invariance j| — > ^ we can deduce the high temper- 
ature behaviour from the low temperature expansion (|4l|) . Using (|35D and 
( flip one gets 

— — = Re ^ -(1 + i— ) + — -^Re * ~(1 + «— ) 

^ 2tt l 2 v tt ;j 48tt 3 L 2 v tt 

' {i -^""[1(1 + ^)] + ... (42) 



11520vr 5 L 2 V 7T 

The expansion in ([12|) is now valid for small (3. For very high temperatures 
(3 « - it is senseless to perform asymptotic expansions of the digamma 
functions in (|42|) for large values of the argument. Rather one could perform 
a Taylor expansion leading to the result 

<9/z 2tt V 2 ; 48vr 3 k 2 j 1 J 768vr 5 V L 15 J 

(43) 

Alternatively one may work with the full expression (|32|). For high tem- 
peratures one may expand the factor factor (irt / k (3) /sinh(7rt / k/3) in powers 

_ 2vrt 

of e ^. One gets 

OA 1 ~ „ T/r l M Ti(2n + 1) . . 
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For small j3 we can use the asymptotic expansions for ty'(z) for large z and 

dA 



thus obtain an expansion for ^ in powers of (3 2 . The result agrees entirely 



with the expansion ([43]) 

The expression for the entropy is obtained by integrating the expression 
and using 



S = + A + + 5- 4 + -1 + ■ • • (45) 

Note that ^"(l) an d ^""(\) are negative so that the leading contribution to 
the entropy is positive. 

One important feature of the high temperature limit is that there is no 
term proportional to the "volume" , which is log \i in this model. 

Contrary to the results of the genus expansion the entropy has a regular 
behaviour at high temperatures. However the specific heat C v = — /3|§ is 
negative. This means that the system is unstable. We should not really trust 
the above thermodynamic expressions for all temperatures. 

We believe that this instability is related to an inherent nonperturbative 
instability of the model defined naively as an inverted harmonic oscillator. 
In fact it is quite unclear how should one define the matrix model so that it 



satisfies all the basic physical requirements [p3J . In a properly defined model 
there should be no instability of the kind discussed above. One possible 
scenario could be : the specific heat should remain positive, but the entropy 
saturates at high temperatures. It is also possible that there is a phase 
transition in the singlet sector itself at (3 ~ -. Such a phase transition is 
distinct from the usual KT transition in this model which is driven by the 
non-singlet states and which is a perturbative phenomenon. 

What is clear from the above discussion, however, is that the genus ex- 
pansion is a bad guide to the behaviour of the thermodynamic entropy at 
high temperatures and that nonperturbative effects are of crucial importance 
in a discussion of the geometric entropy. 

Even if we obtain a physically reasonable expression for the thermody- 
namic entropy it is still not clear how one could use this result to obtain 
the geometric entropy. This is because in this theory interactions are not 
translationally invariant and it is far from obvious how one could extract an 
entropy density which we expect to be position dependent as well. From the 
point of view of the direct calculation of the geometric entropy discussed in 
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the earlier part of this paper this issue is related to the fact that the entropy 
would depend on the region of space which is integrated out to obtain the 
density matrix. 

To obtain the geometric entropy of the underlying string theory, one 
further needs to address the question of the exact correspondence of the 
collective field and the massless scalar of the string theory. The point is, the 
massless scalar of the string theory seems to be related to the collective field 
in a non-local (in space) way |23| . Hence the geometric entropy obtained by 
integrating out the string theory scalar field in some region of the liouville 
space is very different from that obtained by integrating out the fermions 
or the collective fields in some region of the x or r space. Nevertheless our 
results strongly indicate that in all these quantities nonperturbative stronng 
coupling effects play a crucial role. 

Finally to really understand black hole entropy in this string theory one 
needs a clear understanding of the black hole solution in the matrix model. 
There have been several suggestions about this f25j , but the situation is rather 



unclear. In particular there is no contact with known black hole solutions 
of the low energy effective field theory with nontrivial tachyon backgrounds 



2c]. All these issues are intimately tied with a deeper understanding of the 
nature of space of time in the matrix model. 

Acknowledgements: I would like to thank S. Govindarajan, G. Mandal, 
A. Sen, S. Shenker, S. Wadia and especially S. Mathur for discussions. I 
would like to thank A. Jevicki for many discussions and collaboration during 
the initial stages of this work. 



References 

[1] L. Bombelli, R. Koul, J. Lee and R. Sorkin, Phys. Rev. D 34 (1986) 373 

[2] M. Srednicki, Phys. Rev. Lett. 71 (1993) 666. 

[3] C.G. Callan and F. Wilczek, Phys. Lett. B333 (1994) 55. 

[4] D. Kabat and M.J. Strassler, Phys. Lett. B329 (1994) 46. 

[5] J.S. Dowker, University of Manchester Preprint, [hep-th/9401125 . 

[6] F. Larsen and F. Wilczek, [hep-th/9408089 



17 



J.D. Beckenstein, Phys. Rev. D7 (1973) 2333; Phys. Rev. D9 (1974) 
3292 

S.W. Hawking, Phys. Rev. D14 (1976) 2460; Comm. Math. Phys. 43 
(1975) 199; Comm. Math. Phys. 87 (1982) 395; Phys. Rev. D18 (1978) 
1747. 

G. Gibbons and S. Hawking, Phys. Rev. D15 (1977) 2752. 
G. 't Hooft, Nucl. Phys. B256 (1985) 727. 

L. Susskind, Rutgers Preprint, |hep-th/9309145| ; L. Susskind and J. 
Uglum, Phys. Rev. D50 (1994) 2700; D. Lowe, L. Susskind and J. 
Uglum, Phys. Lett. B327 (1994) 226. 

E. Keski-Vakkuri and S. Mathur, Phys. Rev. D50 (1994) 917; T. Fiola, 
J. Preskill, A. Strominger and S. Trivedi, Phys. Rev. D50 (1994) 3987. 

J. Barbon, Phys. Rev. D50 (1994) 2712. 

A. Dabholkar, |hep-th/ 940915^ 

A. Dabholkar, [hep-th / 9408098| ; D. Lowe and A. Strominger, [hep- 



th/9410215 . 



E. Elizalde and S.D. Odintsov, |hep-th/941i024 



J. Atick and E. Wiiten, Nucl. Phys. B310 (1988) 291. 

For reviews and references to the original literature, see I. Klebanov, in 
String Theory and Quantum Gravity, (Trieste Spring School proceed- 
ings, 1991), [hep-th/9108019| ; S.R. Das, in String Theory and Quantum 
Gravity (Lectures in Trieste Spring School, 1992) JEep-th/9211085| ; P. 
Ginsparg and G. Moore, in Recent directions in Particle theory (Pro- 
ceedings of TASI Summer School, Boulder, 1992), |hep-th/93040TT| J. 
Polchinski, frep-th/9411028 



[19] See e.g. K. Huang, Statistical Mechanics, J. Wiley (1963). 
[20] see e.g. S.R.Das in || 



18 



[21] D. Gross and I. Klebanov, Nucl. Phys. B344 (1990) 475; Nucl. 
Phys. B354 (1991) 459. 



[22] M. Bershadsky and I. Klebanov, Phys. Rev. Lett. 65 (1990) 3088. 

[23] J. Polchinski, Nucl. Phys. , B424 (1994) 137; |hep-th/9409i68| ; pej> 



th/941102S . 



[24] M. O'Loughlin, |hep-th/9501039 



[25] E. Martinec and S. Shatashvilli, Nucl. Phys. B 368 (1992) 338; S.R. 
Das, Mod. Phys. Lett. , A8 (1993) 69 and Mod. Phys. Lett. , A8 (1993) 
1331; A. Dhar, G. Mandal and S.R. Wadia, Mod. Phys. Lett. , A7 (1992) 
3703; Mod. Phys. Lett. , A8 (1993) 1701; J.G. Russo, Phys. Lett. B300 
(1993) 336; A. Jevicki and T. Yoneya, Nucl. Phys. , B411 (1994) 64. 

[26] S. Kalyana Rama, Phys. Rev. Lett. , 70 (1993) 3186; Nucl. Phys. B407 
(1993) 191; Phys. Lett. , B321 (1994) 358; A. Peet, L. Susskind and L. 
Thorlacius, Phys. Rev. , D48 (1993) 2415. 



19 



